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Abstract 

Assuming that the universe contains a dark energy fluid with a constant linear equa- 
tion of state and a constant sound speed, we study the prospects of detecting dark energy 
perturbations using CMB data from Planck, cross-correlated with galaxy distribution maps 
from a survey like LSST. We update previous estimates by carrying a full exploration of the 
mock data likelihood for key fiducial models. We find that it will only be possible to exclude 
values of the sound speed very close to zero, while Planck data alone is not powerful enough 
for achieving any detection, even with lensing extraction. We also discuss the issue of initial 
conditions for dark energy perturbations in the radiation and matter epochs, generalizing the 
usual adiabatic conditions to include the sound speed effect. However, for most purposes, 
the existence of attractor solutions renders the perturbation evolution nearly independent of 
these initial conditions. 
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1 Introduction 



The biggest problem in cosmology today is the understanding of the accelerated expansion of the 
universe. Although one could try to attack this question leaving aside the cosmological principle 
or modifying Einstein's gravity, the most classical approach consists of assuming a perturbed 
FLRW universe with a negative pressure component. The minimal model (in terms of number 
of free parameters) compatible with the current data is a cosmological constant, which should 
be perfectly homogeneous by definition. Other candidates (which may or may not alleviate the 
fine-tuning and coincidence problems of the cosmological constant) include, for instance, scalar 
field models, or effective descriptions in terms of a fluid with free parameters yet to be measured. 
A canonically kinetic normalized scalar field would fluctuate, but since in that case the sound 
speed (? s = 5p/8p (computed in the rest frame of the scalar field) is equal to one, local pressure 
would prevent density contrasts to grow significantly. In an effective fluid description, the sound 
speed is a free parameter, and dark energy clustering can be more efficient in the limit in which 
overdensities are not balanced by local pressure perturbations (c s — s> 0) . 

Generally speaking, the study of small perturbations could be used as a tool for discriminat- 
ing between various models with a negative pressure component (cosmological constant, dark 
energy fluid, quintessence or k-essence fields, coupled dark energy, etc.) or a modified theory of 
gravity. One of the major difficulties comes from the fact that the expansion history predicted 
by a given Lagrangian theory of gravity can be reproduced in General Relativity by a dark fluid 
having an appropriate (time varying) equation of state, or by a scalar field with an adequate 
kinetic term and potential. Fortunately, a very precise measurement of clustering properties in 
our universe could at least help to discard some models in favor of others at the level of per- 
turbations [H[2]- However, the spatial fluctuations of typical dark energy models are very much 
suppressed with respect to those of dark matter, and detecting their effect is a real challenge. 

The effects of quintessence perturbations (for which c s = 1) on the CMB and LSS power 
spectra were discussed in [3]. For an (uncoupled) dark energy fluid, there have been several 
studies on the possibility to measure c s , but its value remains unconstrained with present data 
(see e.g. [H EJ [6J [7J El [lOl [TT] ) . In a recent attempt, the authors of [10] used present Cosmic 
Microwave Background (CMB), Large Scale Structure (LSS) and supernovae data (including 
CMB x LSS cross-correlation), and showed that it is possible to see some weak preference for 
c s ^ 1 , but only for a certain kind of early dark energy model in which the equation of state is not 
constant. We may still hope to discriminate between different values of c s using combinations 
of future CMB data with 3-dimensional galaxy clustering data [12] , with CMB x LSS cross- 
correlation data [13] . or with results from a large neutral hydrogen survey such as that of the 
SKA project [14] (see also [i~5]). 

In this work, we focus on an effective description which has already been studied by several 
authors: namely, a dark energy fluid with a linear equation of state p = wp , a constant equation 
of state parameter w close to — 1 and a constant sound speed defined in the range < c 2 s < 1 . In 
Section [2] we review the concept of sound speed for a cosmological fluid. The purpose of Section 
[3] is to clarify the non-trivial issue of initial conditions for dark energy perturbations in the 
radiation era, which are a priori non-adiabatic since (? s > w . We present for the first time the 
initial conditions fulfilled by the dark energy fluid in the synchronous gauge (i.e., in the gauge 
used by most Boltzmann codes), when all other fluids have adiabatic primordial perturbations. 
In Section Owe study analytically the evolution of these perturbations during the matter epoch. 
We derive approximations for the attractor solutions followed by dark energy perturbations 
(both in the Newtonian and synchronous gauges). These new results can be used in the future 
for analytical estimates of the impact of dark energy on structure formation. In Section in 
order to update the analysis of [13], we carry a full Monte-Carlo exploration of the likelihood 
of future mock CMB and LSS data, in order to infer the sensitivity of these data to the dark 



energy sound speed, and to investigate possible parameter degeneracies. Finally, in Section (6] 
we summarize our findings. 



2 The sound speed 

At the level of inhomogeneities, the sound speed of a cosmological fluid plays a similar role to 
that of the equation of state for the background cosmology, and relates the pressure and density 
perturbations as: 

«■' = &■ M 

Defined in this way, the sound speed is gauge dependent. Indeed, the quantity that can be 
assumed to be a definite number (depending on the microscopic properties of the fluid) is the 
ratio 5p/Sp evaluated in the fluid rest frame, often denoted as c s 2 . In another arbitrary frame, 
5p/5p gets corrections related to the velocity of the fluid in that frame, such that in Fourier 
space [I] 

p-Hp = c s 2 5 = c 2 5 + 3H(l + w) {c s 2 - c 2 ) p , (2.2) 

where 5 = 5p/p is the relative density perturbation, 6 the velocity divergence of the fluid, k 
the comoving wavenumber, T~L the conformal Hubble parameter dlna/dr, and c 2 the adiabatic 
sound speed of the fluid. The latter is defined as the (time dependent) proportionality coefficient 
between the time variations of the background pressure and energy density of the fluid, 

p = c 2 p. (2.3) 

For non-interacting fluids one finds that 

w = 3(1 + w) (w - c 2 ) U , (2.4) 

which implies that c 2 = w if w is constant and different from — 1 . In the rest of this work, we 
will focus on the case in which w can indeed be approximated as a constant in time. Besides, we 
will restrict our analysis to < c 2 < 1 . The lower bound prevents dark energy fluctuations from 
growing exponentially, which can lead to unphysical situations; and the upper one is imposed in 
order to avoid superluminal propagation. For a study on the implications of the sign and value 
of (? s for quintessence see [16] . For an expression relating the time evolution of w and c s through 
the intrinsic entropy contribution to the pressure perturbation see [33J (and also 



3 Initial Conditions 

In order to compute the CMB and LSS power spectra, one needs to evolve cosmological per- 
turbations starting from initial conditions deep inside the radiation epoch and far outside the 
Hubble radius. Initial conditions for photons, neutrinos, cold dark matter and baryons are 
reviewed in |17| in the synchronous and Newtonian gauges (sec also [18J). Here, wc want to 
extend this set of relations to dark energy perturbations, especially in the gauge used by most 
Boltzmann codes: namely, the synchronous gauge. Surprinsingly, this issue has been overlooked 
in the literature, without a clear justificatiorH. In practical terms, initial conditions for dark 

3 There were however various studies closely related to this issue. For instance, in [19] various possible initial 
conditions for the quintessence field, and their impact on the CMB were discussed. In [20] ■ the initial conditions 
for a dark energy fluid with quintessence-like perturbations were obtained in a gauge invariant formalism. In [21] , 
the technique of [5D] was extended to interacting dark energy models. In all these works, the dark energy sound 
speed c s was kept fixed to one. 



energy perturbations are essentially irrelevant for most purposes because of the existence of an 
attractor. We wish to clarify this issue and write down explicitly this attractor solution. 

Adiabatic initial conditions on super-Hubble scale derive from the generic assumption that 
for each component i, the density Pi(r, x) can be written as pi(r + 6t(x)) , where Pi(r) stands for 
the background density, and 5r{x) is an initial time-shift function independent of i. Similarly 
the pressure would read pj(r + 5t(x)). Such conditions can be easily justified in all models in 
which primordial perturbations are generated from a single degree of freedom (like the inflaton), 
and/or in cases in which all components have been in thermal equilibrium in the early universe 
with a common temperature and no chemical potential. This form implies 



and hence 5pi/5pi = Pi/p~i , i.e. the sound speed of each species must be equal to its adiabatic 
sound speed. This generic assumption also implies that the total ratio <5p«]/Ei $Pi] ls 
independent of the purely spatial coordinates x . Finally, if all the components do not interact, 
we conclude that 



The fact that for all species the ratios + Wi) are equal to each other is a well-known 

property of adiabatic initial conditions. The meaning of such a relation is not so clear when one 
introduces a dark energy fluid, for which c s 2 > (? a . Hence, in most frames, one has (? s 7^ c 2 . and 
the fluid cannot obey simultaneously (|3.1|) and (|3.2p . This raises the issue of defining sensible 
initial conditions for a dark energy fluid. However, during radiation and matter domination, 
dark energy perturbations tend to fall inside the gravitational potential wells created by the 
dominant component and not much concern has been raised concerning their initial conditions. 
In other words, there is an attractor solution for dark energy perturbations, and their initial 
values are almost irrelevant in practice, provided that for each Fourier mode the attractor is 
reached before dark energy comes to dominate (i.e. provided that initial conditions are imposed 
early enough, and that initial dark energy perturbations are not set to dramatically large values). 
For this reason, in a Boltzmann code like camb [22], initial dark energy perturbations are set 
by default to zero. 

We will derive in the next subsections the attractor solution for a dark fluid with constant 
w and arbitrary c s , assuming that other quantities obey the usual adiabatic initial conditions, 
and in two gauges: the Newtonian and the synchronous ones. For a detailed account of the 
construction of the two gauges and the relations among them we refer the reader to p2] . In 
what follows, we will denote quantities corresponding to the conformal Newtonian gauge with 

(c) 

a superscript. For instance, o x denotes the relative dark energy density in that gauge. No 
superscript will be used for quantities in the synchronous gauge. The transformation equations 
between the two gauges will be summarized below, later in this section. 

3.1 Synchronous gauge 

Early in the radiation era, the total energy density of the universe can be approximated by the 
sum of photon and neutrino densities, with a constant ratio R u = p~ v j (p~ u + p 7 ) . In order to 
find the perturbation evolution on super-Hubble scales and for adiabatic initial conditions, one 
can combine the Einstein, photon and neutrino equations into a fourth order linear differential 
equation for the trace part h of the metric perturbations in Fourier space [T7] . The fastest grow- 
ing mode among the four possible solutions, h ~ {^t) 2 , corresponds to the growing adiabatic 




(3.1) 
(3.2) 



Pi = Pi+pi5r(x) 




(3.3) 



mode. 



Since these conditions are established when the perturbations are still in the super-Hubble 
regime, the product (fcr) <C 1 can be used as an expansion parameter for the solutions of the 
dynamical equations. At leading order 



-\h= -\c{krf = 5 c = 5 b = \b v = , (3.4) 

where the subscripts refer to cold dark matter, baryons, neutrinos and photons; and C is a 
constant. As usual, in order to fully fix the gauge, we impose not only synchronous metric 
perturbations, but also that dark matter particles have a vanishing velocity divergence 8 C . The 
continuity and Euler equations for the dark energy fluid read 

4 = -(i + w) (e + 1 J - 3(c s 2 - w)ns x - 9(1 + w ){c s 2 - c a 2 )n 2 ^ , (3.5) 

c 2 k 2 

9 X = -(l-3c s 2 )H9 x + -^—S x -k 2 a x . (3.6) 

1 + w 

These equations are very general, since the only underlying assumption is that the fluid is non- 
interacting, and allow for the presence of shear stress a x , non-adiabatic sound speed, and a time 
varying w . From now on we assume that the fluid is shear free and has a constant equation of 
state. If the energy density of dark energy at early times is negligible, the solution for the metric 
perturbation h will not change. In order to find the attractor solution, we just need to replace 
h in (13. 5p according to (|3.4j) . and solve equations (|3.5p . (13. 6p . As expected, we find that the 
solution of the homogeneous equation becomes negligible with time, while 5 X and 9 X are driven 
to 

n a — v 2 

2 4 — 6w + 3c s 
a - ;(krfk, (3.8) 



2 4 - 6w + 3c. 



2 ' 



at lowest order in (fcr) . 

This attractor solution does not look like usual adiabatic initial conditions because, in gen- 
eral, c 2 is different from c 2 in the synchronous gauge. Hence, 5p x /5p x cannot be equal to p x /p x . 
However, this solution gives the correct behavior of dark energy perturbations when the other 
components obey adiabatic initial conditions, once the attractor has been reached. Therefore 
they could be called "generalized initial adiabatic conditions". These conditions are valid not 
only for dark energy fluids (w < —1/3) but also for any other fluid with constant w and a = . 
For instance, one can easily check that the usual adiabatic initial conditions for matter and 
radiation can be recovered from (|3.7l I3.8|) by choosing w = c s 2 . 

In a Boltzmann code like camb [22J, the quantitites 5 X and 6 X are fixed to zero at initial 
time for simplicity. For most practical applications, this arbitrary choice does not introduce 
any mistake in the final results, since S x and 6 X are quickly driven to the attractor solutions of 
eqs. (]3.7( 13. 8p . We illustrate this in Figured! for a very large wavelength mode. We suggest 
however to implement eqs. (|3.7l I3.8[) directly into camb's initial condition routine (as we did in 
Section [5] of this work), since this is completely straightforward, and since it offers a guarantee 
that final results are independent of the early time at which initial conditions are defined. 



3.2 Conformal Newtonian Gauge 

The equations that give the evolution of a generic fluid in the conformal Newtonian gauge are 
(see also [TP] ) 



5^ = -(l + w )^) -3<p) -Z{c s 2 - w )m^ -9{l + w){c s 2 -Ca 2 )^ 2 -^, (3.9) 

0(«O = -(l-3c s 2 )nd^ + ^-S^-k 2 a^+k 2 i;. (3.10) 
The metric perturbations 

(j) = n-aU, (3.11) 

tj) = a + aH, (3.12) 

can be obtained from those of the synchronous gauge using 

2/c 2 a = A + 6r) , (3.13) 

where 77 is the traceless part of the metric scalar perturbation in the synchronous gauge in 
Fourier space. Using these last equations one can immediately check that the product aH has 
zero time derivative, so that 4> and ip are time independent at lowest order in (kr) . One can 
either solve directly (13.91 13.10]) . or use our results <\3.7\ I3.8j) for the behavior of S x and 6 X in the 
synchronous gauge, and perform the gauge transformations 

4 C ) = -3(1 + w)aH + S x , (3.14) 
0( c ) = ak 2 + 9 x , (3.15) 
4 C) = <*x, (3.16) 

that are valid for non-interacting fluids. The two methods give 

4 C) = ~(l + w)^ + S x , (3.17) 
6$ = ^(kr)k + 9 x , (3.18) 

where 

90 

1P = C (3.19) 

and 

Notice that the leading contributions to the velocity and density perturbations in the conformal 
Newtonian gauge are independent of the speed of sound and such that 

( c ) 

5x - _!L - AW - A( c ) - 3 r(c) _ K(c) ( n 2l) 

l + w ~ c ~ b ~4 U ' 1 ' 

Therefore, the usual adiabatic conditions are recovered, and c 2 enters only in the next order 
correction^. Indeed, the Newtonian gauge is the one in which, beyond the Hubble scale, c 2 



4 Note that the gauge transformation law (|3,14[) implies that 



1 + W 1 + Wj 1 + W 1 + Wj 



(3.22) 



is equal to c 2 at leading order, even when c s 2 ^ c 2 . This can be checked by keeping the 
dominant terms in (|3.14l 13. 15H . and replacing 6 X and X by these values in (12. 2|) : one gets 
c 2 6 x = c 2 6 x c ^ + O(kr) 2 . So, in the Newtonian gauge and on super-horizon scale, 5p/5p is equal 
to p/p for any fluid with contant w, and the common intuition according to which super-Hubble 
fluctuations behave in the same way as background quantitites is recovered. 



4 Late time attractors 



In the previous section, we found the attractor solutions for dark energy perturbations during 
radiation domination. Here we will derive similar solutions during matter domination. These 
results can be used to provide initial conditions for dark energy perturbations in problems in 
which following the behavior of cosmological perturbations during radiation domination is not 
relevant. 

If we assume that the energy density of photons and (massless) neutrinos is negligible deep 
inside matter domination, the perturbations can be studied using a two-fluid approximation. 
One of the fluids is formed by baryons and cold dark matter (which cannot be distinguished 
from each other) and the other one is dark energy. This description in terms of two components 
can be accurately used to study the growth of matter perturbations up to the present day. 
Mathematically, the problem consists of a system of six independent equations with six variables: 
the density and velocity perturbations of the two fluids plus the scalar metric perturbations. We 
will now find the relevant growing modes of the perturbations in the two gauges and conclude 
this section with some remarks concerning the initial conditions for dark energy. 



4.1 Synchronous gauge 

In the synchronous gauge, since we consider that cold dark matter and dark energy are shear 
free, the Einstein equations imply that the two metric degrees of freedom rj and h are related to 
each other. Eliminating r/ in terms of h, and expressing h in terms of 5 C (the continuity equation 
for cold dark matter gives h = —25 c ) , we obtain a system of two reasonably short second order 
differential equations that describe the evolution of density fluctuations [23] : 



6 C + H5 c --U 2 n c 5 c 



[l + 3c s 2 )6 x + 9(l + w)H(c s 



w 



k 2 



(4.1) 



+ [3 



w 



)H-T(k,H)]5 x 



w 



[1 + 3wQ x - 6c s 2 ) H 2 + ZF(k, U)U 



w 



{l + w)6 c -{l + w)JF{k,H)5 c , 



(4.2) 



for any fluid j with constant equation of state Wj. Since the usual adiabatic conditions hold at leading order 
in the Newtonian gauge, one may naively infer from the above equality that they hold also in the synchronous 
gauge. This is not correct since on the left-hand side, the two terms are dominated by order zero terms in a (kr) 
expansion, while on the right-hand side the leading terms are of order two. Assuming that the fluid j has an 
adiabatic sound speed (like cold dark matter or photons), a full order-two calculation of all the terms leads to 



t (c) 



5f 



1 + W 1 + TO,' 1 + W 1 + Wj 



c 



3(3 



4 - 6w + 3cj 



(kr) 



(3.23) 



The righ-hand side does not vanish since (c 2 a — w) is by assumption strictly positive. Being of order two, this 
difference contributes to the solutions at leading order in the synchronous gauge, but only at next-to-leading order 
in the Newtonian gauge. 



where 



and 



(1 + w) % = vKW) ("^ + (1 + w)Sc ~ 3 ^ ~ ^ H5x ) ' (43) 



V(k,U) = k 2 + 9(c s 2 - w)U 2 , (4.4) 



T{k,U) = -9(l + 3wn x )^^H 3 -(l-3c s 2 )H. (4.5) 

In the purely matter dominated epoch (f2 c = 1, Q x = 0) the equation that describes the 
evolution of matter perturbations is the classical growth formula: 

s c + ns c - ^n 2 n c 5 c = o . (4.6) 

Its general solution is a linear combination of a growing mode (S c ~ a) and a decaying one 
(<5 C ~ a -3 / 2 ) . In order to find the relevant attractor solution for dark energy perturbations, one 
should keep the first of these two solutions 

5 C oc (kr) 2 . (4.7) 

The behavior of cosmological perturbations depends on their wavelength as compared to the 
characteristic scales of the problem. In our case, and in terms of comoving scales, there are two 
relevant quantitites: the comoving Hubble scale 'H -1 (which gives the order of magnitude of 
the causal horizon associated to any process starting after inflation), and the comoving sound 
horizon %J l = c s 1-L~ l . Analytical approximations for the evolution of dark energy perturbations 
can be obtained in the three regimes defined by these two scales. 

For super-Hubble perturbations with k <C TL , one can approximate (|4.5I) as J-(k, 7i) ~ 
(3c s 2 — 2) H , and the homogeneous part of (|4.2p gives two decaying modes for 5 X . The growing 
solution, sourced by the dark matter perturbations, is 



5 - 6c s 2 
15w + 9c, 



5 x = (l + w)- ~~~^~t~2^£ ' k^U. (4.8) 



If instead T~L S S> k S> H , the perturbations are below the Hubble scale but above the sound 
horizon, and J-(k, Ti) ~ (3c s 2 — l) % . As in the previous case, the growing mode for the density 
perturbations of dark energy comes from the one of dark matter: 

5 x = (l + w) l ~ 2Cs 2 S C , % s ->k->%. (4.9) 
1 — 3w + c s 

These formulas agree with the results of [23] (when transformed into the synchronous gauge) 
if we take the sound speed to be zero. The equations (|4,8p and (|4.9p can be used to define 
initial conditions for a dark energy fluid during matter domination, provided that at the initial 
time all relevant wavelengths are still larger than the sound horizon. Below the sound horizon, 
it becomes more difficult to obtain simple analytic approximations. When T~L S <C k , we can 
still approximate F(k,'H) as in the previous case, but the term proportional to c s 2 A; 2 in (|4,2p 
becomes dominant. One gets: 

5 X + (1 - 3w)H5 x + c 2 k 2 b x = - (1 + w)(l - 2c s 2 )U 2 6 c , U s < k . (4.10) 

The general solution of the homogeneous equation goes as t?>w-i/2 w [fa a multiplicative factor 
that is a linear combination of the Bessel functions J Sw _ i/ 2 (c s kr) and Y 5w _ 1 / 2 (c s kT) . The 



particular solution for the complete equation including 5 C can also be written in terms of non- 
elementary functions. Dark energy perturbations are anyway suppressed with respect to dark 
matter ones in this regime, since below the sound horizon (i.e below the Jeans length of the fluid) 
the pressure perturbation can resist the gravitational infall. Indeed, equation (|4.10j) implies that 
in the limit k — > oo , the dark energy density contrast 6 X should vanish. 

4.2 Conformal Newtonian gauge 

In the conformal Newtonian gauge, the complete equations for the evolution of density pertur- 
bations, equivalent to (|4.ip and (|4.2p . become longer, because metric perturbations cannot be 
trivially replaced in terms of 5 C . Having obtained the solutions in the synchronous gauge, it is 
simpler to apply (|3.14p rather than solving the conformal Newtonian equations directly. In the 
limit of pure matter domination, i.e. H x = , one can easily prove that 

h = 2k 2 a, (4.11) 

and therefore 

= (l + 3^)* c , 
4 C > = 5 X + 3(1 + w)^5 c 
These equations show that outside the Hubble radius, the solution is driven as usual to 

4 C > = (1 + «/)#>. (4.14) 

The same equations also imply that in the regime 7i s ^> k S> H , (|4.9p remains valid in the 
conformal Newtonian gauge (relations valid inside the Hubble radius are expected to be gauge 
independent). 

4.3 Attractor solutions 

An important feature in the evolution of dark energy perturbations in the two gauges, which 
is common to the three regions we have studied, is that initial conditions for the dark energy 
perturbations are almost irrelevant for the evolution in the purely matter dominated epoch. 
During this period, dark energy fluctuations track matter inhomogeneities. Above the sound 
horizon, as soon as the attractor solution of (|4.8|) or (|4.9p is reached, the ratio S x /S c only depends 
on the sound speed and equation of state of the dark energy fluid. Notice that in the newtonian 
gauge the sound speed dependence actually disappears outside the Hubble scale. 

In the dark energy dominated period (£l x — > 1 ) , the roles of the fluctuations are exchanged 
and matter perturbations are sourced by dark energy ones. However, since we have seen that 
it is only the initial value of 5 C that determines 5 X in the matter epoch, the initial conditions 
for the evolution in the dark energy period are in reality given only by 5 C . In accordance with 
the results of section [3], this argument can be extended back into the radiation era. Besides, 
from the previous reasonings, it is clear that the velocity perturbations are also unimportant. In 
conclusion, the amount of dark energy perturbations today can be well estimated just by knowing 
the dark matter perturbations at some initial time in the radiation or the matter epoch. We 
illustrate this behavior in Figured] for a very large wavelength mode. 
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Figure 1: The three thick lines show the evolution of the ratio S x /Sy in the synchronous gauge, 
obtained numerically with camb, in a model with: w = —0.9, (? s = 0.1, standard values of 
the other cosmological parameters, and adiabatic initial conditions for photons, neutrinos, cdm, 
baryons. We choose the case of a very long wavelength mode (k = 2.3 x 10 _6 Mpc _1 ) which 
remains outside the Hubble radius during all relevant stages: radiation domination (RD), matter 
domination (MD) and dark energy domination (DED). We integrated this mode starting either: 

(i) from the "usual initial condition" 5 X /Sj = (1 + w)/(l + io 7 ) with w 7 = 1/3 (upper thin 
horizonal line), which has no physical justification in the synchronous gauge in this context; 

(ii) from eq. A3 . Tj) (middle thin horizontal line) and (I3.8P : (Hi) from 5 X = 0, like in the public 
version of camb. In each case, the solution quickly evolves in such way to fulfill eq. (j3.7|) (middle 
thin horizontal line) during radiation domination, and then eq. (14. 8p (lower thin horizontal line) 
during matter domination. 




5 Detectability 

The detectability of the sound speed of dark energy has already been studied by various authors, 
under different assumptions and for various datasets. For example, for a model with constant w 
and c s (identical to the one we consider in this work), the authors of [TO] showed that the combi- 
nation of present CMB, LSS and supernovae data is not sensitive at all to the dark energy sound 
speed. Since dark energy perturbations would change the growth rate of matter inhomogeneities 
on intermediate scales (between the Hubble radius and the dark energy sound horizon), they 
can affect the CMB photon temperature through the Late Integrated Sachs- Wolfe (LISW) effect. 
Small variations in the LISW effect are difficult to detect in the CMB temperature spectrum, 
due to cosmic variance and to the fact that we only see the sum of primary anisotropics and 
LISW corrections. However, the LISW contribution can be separated from primary anisotropics 
by computing the statistical correlation between CMB and projected LSS maps of the sky. In 
|10j . most of the presently available cross-correlation data were included in the analysis, but 
current statistical and systematic error bars are far too large for probing sub-dominant dark 
energy clustering effects. 

One could think of constraining the dark energy sound speed in the future, either by mea- 
suring with better accuracy the auto-correlation function of matter distribution tracers (galaxy 



surveys, cluster surveys, lensing surveys, . . . ) or, again, by studying their cross-correlation with 
CMB maps in order to extract the LISW contribution. The first option was studied for instance 
in [12], and the second in [13]. In the latter, the authors focus on the cross-correlation of a 
CMB experiment similar to Planck with a survey comparable to the Large Synoptic Survey 
Telescope (LSST) project. The author of [12] considered the combination of Planck with fu- 
ture galaxy redshift surveys, not including any cross-correlation information, but using the full 
three-dimensional power spectrum of galaxies instead of their angular power spectrum sampled 
in a few redshift bins (the loss of information on the LISW effect can then be compensated 
by more statistics on the auto-correlation function). Both works reached the same qualitative 
conclusion that the next generation of LSS surveys could discriminate between quintessence-like 
models with c s = 1 and alternative models with a sound speed c s < 1, and estimated under 
which threshold in c s the discrimination would be significant. 

The works of p2] and [13] are based on analytic estimates of the measurement error for each 
cosmological parameter, using the Fisher matrix of each future data set. This is the second-order 
Taylor expansion with respect to the cosmological parameters of the data likelihood around its 
maximum, i.e. in the vicinity of a putative fiducial model. This method should be taken with a 
grain of salt whenever the dependence of the observable spectrum on a given parameter cannot 
be approximated at the linear level within the range over which this parameter gives good fits 
to the data. This is precisely the case for the model at hand. The dark energy sound speed will 
not be pinned down with great accuracy using the Fisher matrix approach. Its variation within 
the error bar has a non-trivial effect on the total matter power spectrum which, although small, 
cannot be captured at the linear level, since it amplifies perturbations over a range of scales 
depending on c s itself. This means that Fisher matrix estimates for the error Ac, should only 
be considered as a first-order approximation. This caveat was already emphasized in Section 
IV of |13j . Hence, it is worth checking the results of [13] with a full exploration of the likelihood 
of mock Planck + LSST data. We generated such data for some fiducial models including a 
dark energy fluid, sampled the data likelihood with a Monte Carlo Markhov Chain (MCMC) 
approach, and inferred the marginalized posterior probability distribution of the dark energy 
sound speed. Unlike the Fisher matrix approach, this method takes into account the precise 
effect of each parameter throughout the parameter space of the model. Hence, it can assess 
realistically whether the data can resolve non-trivial parameter degeneracies. As in [12], we 
included the total neutrino mass in the list of parameters to measure, in order to check whether 
some confusion between the effect of dark energy clustering and massive neutrino free-streaming 
could arise when the dark energy sound horizon is comparable to the neutrino free-streaming 
scale. Our mock data sets were generated and fitted with two modified versions of the public 
CosmoMC package [26] which were already developed and described by the authors of [27] 
and [28] • 

5.1 Planck with lensing extraction 

In a first step, we focus on the potential of the Planck experiment alone, assuming Blue Book 
sensitivities for the most relevant frequency channels, and a sky coverage of 65% (as summarized 
in [29]). Given that small variations in the LISW effect cannot be probed due to cosmic variance, 
the only hope to detect dark energy clustering with Planck data alone is through lensing extrac- 
tion. A lensing estimator can in principle measure the deflection of CMB anisotropies caused 
by gravitational lenses (e.g. galaxy clusters). This techniques allows to reconstruct the power 
spectrum of the gravitational potential projected along the line of sight, and to disentangle a 
fraction of the total LISW contribution to temperature anisotropies (by cross-correlating the 
lensing map with the temperature map). Like in [27j, we ran CosmoMC in order to fit some 
mock Planck data, including reconstructed lensing data with a noise level corresponding to the 



"minimum variance quadratic estimator" of [30J. The mock data generator and all modifica- 
tions to the CosmoMC package are publicly availabl^E Our results show that Planck alone is 
completely insensitive to the dark energy sound speed for any constant and reasonable value of 
w. Indeed, when assuming a flat prior on log 10 [c s ] in the range [—3,0], we obtained a nearly 
constant posterior distribution 'P(log 1 g[c s ]). We checked this results with several assumptions 
on the fiducial values of w (bewteen —0.8 and 1), c s (between 10~ 2 and 1) and £m„ (between 
0.05 eV and 0.2 eV). This negative result can be attributed, on the one hand, to the limited 
precision of lensing extraction with Planck, and on the other hand, to the fact that the projected 
gravitational potential probed by CMB lensing gives more weight to high redshifts (of the order 
of z ~ 3) than to the redshifts at which dark energy comes to dominate and eventually affects 
the growth of total matter perturbations (z < 1 for models with a constant w). In the rest of 
this section, we will not include Planck lensing extraction data anymore. 

5.2 Planck plus LSST galaxy data 

In a second step, we performed a combined analysis of mock Planck and LSST data, including 
information on the angular auto-correlation function of LSST galaxy maps (in 65% of the sky 
and in 6 redshift bins spanning the range from z = to z ~ 5), and on the cross-correlation 
between these maps and CMB temperature maps. Our assumptions concerning LSST data 
binning, selection function, bias, density and coverage are identical to those in |31|. I28j . The 
corresponding expression for the likelihood of correlated CMB and galaxy data is given in [28] . 
We generated two mock data sets, assuming two fiducial cosmological models with: standard 
values of the six ACDM parameters, a total neutrino mass ^ m v = 0.2 eV (equally split between 
the three species for simplicity), a constant dark energy equation of state parameter w = —0.9, 
and a dark energy sound speed equal to either c s = 1 or c s = 10~ 2 . For each model we ran 
sixteen chains with flat priors on the usual CosmoMC parameter basis (u>b, uj c , 6, r, logL4 s ], 
n s , f u , w) plus log 10 [c s ], imposing a restriction log 10 [c<j] < 0. After reaching a convergence 
criterion R — 1 < 0.01, we obtained the marginalized posterior probability of each of these 
parameters, and of the derived parameters ^2m u (total neutrino mass) and Hq. We show in 
Table Q] the 68% standard error for each of these. In the first plot of Figure [21 we display the 
marginalized distribution of log 10 [c<j] in the two models. The correlation between log 10 [cs] and 
other parameters appears to be small, except for w, f u (or ^2m u ), and 9 (or Hq). In Figure [21 
we also show the two-dimensional 68% and 95% confidence levels contours for these pairs of 
parameters. 

We learn from these runs that correlated CMB+galaxy data from Planck+LSST may dis- 
criminate between various sound speed values with a standard deviation Alog 10 [c s ] of order 1 . 
In other words, such a data combination can give an indication on the order of magnitude of c s . 
With a fiducial value c s fiducial — i ; the lower bound of the 68% (resp. 95%) confidence interval 
is log 10 [c s ] > — 1.1 (resp. — 1.8) . With a fiducial value c s fiducial = 10 -2 , these numbers become 
log 10 [c,s] > —1.9 (resp. —2.3). The 68% confidence intervals always include log 10 [c s ] = 0, al- 
though the case of c s fiducial = 10~ 2 appears to be at the threshold below which the data would 
impose a negative 68% upper bound on log 10 [c s ]. These estimates are not as optimistic as those 
in [13], most probably because of our more accurate technique for error forecast (based on a 
marginalization of the actual likelihood over the full parameter space). They are also less dra- 
matic than those of [12], inferred from the Fisher matrix of a different dataset. Still, they show 
in a robust way that with Planck+LSST one can exclude dark energy models with maximum 
clustering (i.e. c s — > 0), which is not the case with current data (and will not be the case even in 
one year from now with Planck data). This limitation is due to the smallness of the effect and not 
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Figure 2: (Top left) Marginalized posterior probability distribution of log 10 [c s ] for Planck+LSST 
and the two fiducial models described in the text and in Table in which c s = 1 (red lines) 
or c s = 10~ 2 (black lines). The dotted lines show the mean likelihood for comparison. (Other 
plots) For the same models and mock data, joint 68% and 95% confidence contours for log 10 [c s ] 
and the three parameters most correlated to the dark energy sound speed. 



to parameter degeneracies, since logic's] is not particularly correlated with other parameters. 
Of course, the higher is w, the longer will be the dark energy domination and the stronger will 
be the effect of c s . This explains the correlation between log 10 [c s ] = and w seen in Figure [2j 
if w were found to be closer to —0.8 than to —1, the measurement of c s would be considerably 
easier. This is in agreement with the c s dependence of the bounds on w found in [25]. The 
joint (log 10 [c s ], w) likelihood contours actually suggest that even with u; fiducial = -0.85 , a value 
like c s = 10 -2 could be pinned down with a fairly good accuracy. Like in [12] . we find that the 
degeneracy between the dark energy sound speed and the neutrino mass is not significant: this 
can be understood from Figure 5. of [12] . which shows that dark energy clustering and neutrino 
free-streaming can produce a step in the matter power spectrum on similar scales, but with 
different shape. The step induced by dark energy is sharper because this effect occurs during a 
brief period of time compared to non-relativistic neutrino free-streaming. 



6 Conclusions 

In sections [3] and [4J we have studied analytically the behavior of dark energy perturbations 
in the radiation and matter epochs in the synchronous and conformal Newtonian gauges. We 
have written down the formulas that generalize the usual adiabatic initial conditions for small 
inhomogeneities when a dark energy fluid with constant equation of state and speed of sound 
is considered. Our equations (|3.7l \3.8\i can be readily implemented in the initial conditions of 
e.g. camb, which applies to early radiation domination. Instead, equations (I4.8| I4.9D - or their 
counterpart ()4.144 I4.9|) in the Newtonian gauge - provide correct initial conditions above the 
sound horizon for codes simulating only the growth of matter fluctuations during radiation and 
dark energy domination. The fact that there is an attractor behavior for the evolution of the 
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Table 1: Standard deviation on each cosmological parameter for correlated CMB+galaxy data 
from Planck+LSST for two fiducial models with standard values of the six ACDM parameters, 
a total neutrino mass ^2m u = 0.2 eV (equally split between the three species), a constant dark 
energy equation of state w = —0.9 and a dark energy sound speed equal to either c s = 1 or 
c s = 10~ 2 . The last two lines show derived parameters. For log 10 [c s ], instead of quoting the 
standard error, we provide the lower bound of the 68% confidence interval (the upper bound 
being imposed by the prior, log 10 [c s ] < 0). 

perturbations guarantees that dark energy fluctuations are driven to solutions dictated by the 
dark matter ones. This implies that the final state of dark energy today is independent on 
its possible initial conditions, provided that the attractor solution is reached way before dark 
energy comes to dominate the background expansion; using the above set of initial conditions 
guarantees that such a condition is fullfilled. 

After clarifying these issues, we have studied in section [5] the prospects for detecting dark 
energy perturbations using CMB data from the Planck satellite, cross-correlated with galaxy 
distribution maps obtained with an LSST-like instrument. We have chosen to focus on models 
of dark energy consisting in a fluid with constant equation of state w and constant sound speed 
c s . We have performed a full exploration of the likelihood of our mock data using a Monte Carlo 
approach, in order to cross-check previous Fisher matrix estimates in a more robust way. We 
find that the confidence interval for c s inferred from the data will potentially allow us to put 
a lower bound on c s , and to constrain at least its order of magnitude (although with a slightly 
smaller significance than in earlier analytic estimates). Our results build upon previous works 
on the same topic [H [131 fl2| [TO] and complement other researches were different models have 
been studied (see for instance |10j for early dark energy and |32| for the case in which there is 
a coupling to dark matter). 
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